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Abstract. We develop a new method of solving Bethe-Salpeter (BS) equation in Minkowski space. It is 
based on projecting the BS equation on the light-front (LF) plane and on the Nakanishi integral rep- 
resentation of the BS amplitude. This method is valid for any kernel given by the irreducible Feynman 
graphs. For massless ladder exchange, our approach reproduces analytically the Wick-Cutkosky equation. 
For massive ladder exchange, the numerical results coincide with the ones obtained by Wick rotation. 
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1 Introduction 

BS equation £Q is an important tool to study the relativis- 
tic bound state problem in a field theory framework (see 
for review 0). For a bound state of total momentum p 
and in case of equal mass particles, it reads 

i 2 

^ k,P) = [(§ + k) 2 - m 2 + ie] [(§ - kf -m 2 + ie] 

x / (1) 

where <P is the BS amplitude, iK the interaction kernel, 
m the mass of the constituents and k their relative mo- 
mentum. We will denote by M — ^pp 2 the total mass of 
the bound state, and by B = 2m — M its binding energy. 

It was recognized from the very beginning that, when 
formulated in Minkowski space, the BS equation has sin- 
gularities which make difficult to find its solution. These 
singularities are due to the free propagators of the con- 
stituent particles 

G^\k, P ) = G^G^^ 



(| + k) 2 — m? + ie (| — k) 2 — m 2 + ie 

but can also result from the interaction kernel itself. 

To overcome this difficulty, Wick ^ formulated the BS 
equation in the Euclidean space, by rotating the relative 
energy in the complex plane feo — > iko. This "Wick ro- 
tation" led to a well defined integral equation which can 
be solved by standard methods. Most of practical appli- 
cations of the BS equation have been achieved using this 
technique |2j and recent developments make its solution 
a trivial numerical task 0]. Another method - the varia- 
tional approach in the configuration Euclidean space - was 



recently developed in [5]. Whereas the total mass of the 
system is unchanged by the Wick rotation, the original BS 
amplitude is however lost and the "rotated" one can no 
longer be used in calculating other physical observables, 
like for instance form factors. 

Thus, fifty years after its formulation, obtaining the BS 
solutions in the Minkowski space is still a field of active re- 
search. A successful attempt was presented in 6 , based on 
the Nakanishi integral representation of the BS function 
0. However, formal developments displayed in are a 
matter of art and the obtained equation has been derived 
and solved only for the ladder kernel. Another approach 
in Minkowski space for separable interactions was devel- 
oped in |S] and applied to the nucleon-nucleon system. On 
another hand, an equation obtained by projecting the orig- 
inal BS equation on the LF plane, was derived and solved 
in j9]. An approximate LF kernel was there obtained as 
an expansion of the BS one but the original BS amplitude 
has not been reconstructed from its LF projection. 

The aim of this paper and the forthcoming one is to 
present a new method of solving the BS equation without 
using the Wick rotation. Our method is based on an inte- 
gral transform of the initial equation which removes the 
singularities of the BS amplitude. This integral transform 
consists in projecting the BS equation on the LF plane, 
defined by u> ■ x — with lo 2 — ^UJ- The particular 
choice u) = (uiq,u>) = (1,0,0,-1) results in the standard 
LF form t + z = and in the LF projection used in [S]. In 
our approach, the BS amplitude maps onto the LF wave 
function while the transformed equation - in contrast to 
PI - is derived without any approximation. This equation 
remains equivalent to the original BS one, therefore pro- 
viding the same binding energies, and the initial BS am- 
plitude is easily reconstructed from its solution. Although 
results presented here concern only the ladder kernel, our 
method is not restricted to a particular interaction. For 
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more complicated kernels, e.g. the cross box, calculations 
become more lengthy, but the additional difficulties are 
due to evaluating the Feynman diagram itself and not to 
the solution of the equation. 

In order to present the method more distinctly, we con- 
sider the case of zero total angular momentum and spinless 
particles. 

The plan of the paper is the following. In sect. |2 we 
give the integral transform used to project the BS equa- 
tion on the LF plane and we derive a new and equivalent 
equation. In sect. [21 the corresponding ladder exchange 
kernel is calculated analytically. In sect. the numeri- 
cal solutions for the ladder case are found and compared 
to the results obtained using other methods in Euclidean 
space. Sect. contains concluding remarks. Details of the 
calculations are given in appendices El an d 10 The re- 
sults concerning cross box kernels are presented in the 
next paper 



2 Projecting the BS equation on the LF plane 

Our method is inspired by an existing relation between the 
BS amplitude <P(k,p) and the two-body LF wave function 
ip(kx, x). This wave function can be obtained by project- 
ing the BS amplitude on the LF plane. We will apply 
below the LF projection to the BS equation in Minkowski 
space. Though this projection can be considered as a for- 
mal transform, we will start by reviewing its derivation, 
in order to show more clearly how the singular behaviour 
of <P{k,p) gives rise to a non-singular tp(k±,x). 

BS amplitude is defined as the matrix element between 
the vacuum (0| and a state \p) of the time ordered product 
of two Heisenberg operators: 



${x u x 2 ,p) = (0\T {(p(x 1 )ip(x 2 )} \p). 



(3) 



In general, the state vector \p) can be taken in different 
representations. In the LF quantization, it has the form: 
(see e.g. eq. (3.1) from [TH]): 



\p) = J i>(ki, k 2 ,p, wt)2(w • p)5(ki + k 2 — p — lut)cIt 
x (2^ 2 rf3fcl d3k2 a t a t | \ 



+ 



(4) 



where at is the creation operator and = V m 2 + k 2 . The 
two-body Fock component ip is shown explicitly, whereas 
the higher ones are implied. All the four-momenta are on 
the corresponding mass shells kf = m 2 , p 2 = M 2 , (lot) 2 = 
and fulfill the conservation law 

ki + k 2 = p + CJT. 

Projecting the BS amplitude <&{xi,x 2 ,p) on LF plane 
means that its arguments are constrained to u> ■ x\ = ui ■ 
x 2 — 0. Coming to the momentum space, we still keep this 



constrain. Let us evaluate the quantity: 

J(k x ,k 2 ,p) = J d^x^x^uj ■ Xl ) 5{lo ■ x 2 )e l{kl - Xl+k2 - X2] 



(5) 



We substitute here the right-hand side of © with \p) given 
by On the LF plane the Heisenberg field ip(x) in © 
turns into the Schrodinger (free) one, represented as 



<p{x) 



a{e~ lkx 



a k e 



d 3 k 



Then the two-body component ip only survives in \p) and 
J(ki, k 2 ,p) is expressed through it. 

Now express <P(xi,x 2 ,p) in J^J through its Fourier trans- 
form. Translational invariance imposes <P to have the form 



<P(xx,x 2 ,p) 



(2tt) 3 / 2 



<P(x,p) e 



-ip-(xi+x 2 )/2 



where <P(x,p) is the reduced amplitude and x = X\ — x 2 . It 
is expressed through the momentum space BS amplitude: 



d A k 



$(k,p) e 



-ik-x 



where $(k,p) satisfies the BS equation QJ. Substituting 
theses formulas in J5j), we find that J(ki, k 2 ,p) is expressed 
through the integral <P{k + (3uj,p)d(3. Comparing two 
expressions, we obtain the relation |10j : 



i/j(k ± ,x) = 



(u ■ ki){w ■ k 2 ) 
tt(lo ■ p) 



<P(k + (3uj,p)dp. (6) 



In the standard LF approach the /3-integration in © turns 
into the fc_ -integration with fc_ = fco — k z . Wave function 
ip(k±,x) in 10 is parametrized in terms of the standard 
LF variables k± 7 x (see [TU]): 



5lJ_ - xp±, 



LO ■ ki 

LU ■ p 



(7) 



The _L-components are orthogonal to oj. 

LF wave function ip, as any wave function, has no sin- 
gularities in physical domain. Equation JBJl can be thus 
viewed as an integral transformation of the BS amplitude 
leading to a non-singular function. It suggests to apply 
this transformation to the BS equation itself: 



dp<P{k + (3uj,p) = 



d(3G ( Q 12) (k + pLj,p) 



d A k' 



(8) 

■iK(k + Puj,k',p)<P(k',p), 



in order to obtain an equivalent equation free of singular- 
ities. This constitutes the key point of this work. 

Apart from the trivial kinematical factor ( CJ,fc i)("' fc2 ) 
the left-hand side of is the LF wave function -0, eq. 
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10, whereas, the right-hand side still contains the "non- 
projected" BS amplitude 4>(k',p). To make the LF wave 
lunction appear explicitly in the right-hand side too, and 
thus formulate an equation in terms of tp, would need to 
invert equation JfJJl. Instead, we substitute, in both sides of 
the equation (JSJ), the BS amplitude in terms of the Nakan- 
ishi integral representation 0Cj : 



<P(k,p) = 



— , dz' / drf' 

V47T 7-1 Jo 



(9) 



[7' + m 2 - \M 2 - k 2 - p ■ k z' - ie] 



3 ■ 



In more general form of this representation, the denomi- 
nator appears in the degree 2 + n, where n is a dummy 
integer parameter. For simplicity, we chose here its mini- 
mal value n — 1 . Greater value of n may result in a more 
smooth solution 6 . 

A similar representation exists for non-zero angular 
momentum. It is valid for rather wide class of the solutions 
which are consistent with the perturbation-theoretical an- 
alyticity. This leads (see appendix^for the detail of calcu- 
lations) to the following equation for the weight function 
9{l, 2): 



7' + 7 + z 2 m 2 + (1 — z 2 )k 2 



dry' / dz 1 V{"i,z;i,z')g{i,z'), 
J— 1 



(10) 



This is just the eigenvalue equation of our method. It is 
equivalent to the initial BS equation JTJ- The total mass 
M of the system appears on both sides of equation ltTU|l 
and is contained in the parameter 



-M 2 



(11) 



As calculations jS] show, 5(7, z) may be zero in an interval 
< 7 < 7o- The exact value where it differs from zero is 
determined by the equation (|10|l itself. 

The kernel V, appearing in the right-hand side of eq. 
(|10|) . is related to the kernel iK from the BS equation by 



1/(7, z; 7', z) 



u> ■ p 



-il{k + (3uj,p)d(3 



[(§ + k + f3uj) 2 - m 2 + ie] 



[(§ - k - /3lu) 2 - m 2 + ie] 



with 



I(k,p) = J 



d 4 k' 



iK(k,k',p) 



( 27r ) 4 [fc' 2 + p ■ k'z' 



(12) 



(13) 



The particular cases of the ladder kernel and of the Wick- 
Cutkosky model |3lll2| are detailed in the next section and 
for the cross ladder kernel - in the next paper Once 
3(7, z) is known, the BS amplitude can be restored by eq. 
©• 

The variables (7, z) are related to the standard LF 
variables (JZJ)as7 = fc 2 _,2: = l — 2x. The LF wave function 
can be easily obtained by 



ip(k_L,x) 



I 



x(l - x)g(y,l- 2x)di 



7' + k\ + m 2 - x(l - x)M< 



(14) 

Eq. I|10|) can be transformed to the equation for the LF 
wave function ip(k±,x) (eq. (|32|l in appendix 0), though 
this requires inverting the kernel in the left-hand side of 
The initial BS equation JHJ, projected on the LF 
plane, can be also approximately transformed (see ap- 
pendix 0) to the LF equation: 



kl 



M ip{kx,x) 



x(l — x) 

^ J ^(k' ± ,x')V LF (k' ± ,x';k ± ,x,M 2 



(15) 



d 2 k' ± dx' 
2x'(l - x') 



with the LF kernel Vlf given, for ladder exchange, by eq. 
(I5*5|l in appendix 151 

It is worth noticing that the LF wave function (|14|l 
is different from the one obtained by solving the ladder 
LF equation 1)15(1 . as it was done e.g. in ref. ^21- The 
physical reason lies in the fact that the iterations of the 
ladder BS kernel (Feynman graph) and the ladder LF ker- 
nel (time-ordered graphs) generate different intermediate 
states. The LF kernel and its iterations contain in the 
intermediate state only one exchanged particle, whereas 
the iterations of the ladder Feynman kernel contain also, 
many-body states with increasing number of exchanged 
particles (stretched boxes). This leads to a difference in 
the binding energies, which is however small P~3;. For- 
mally, this difference arises because of the approximations 
- explained in appendix 151 - which are made in deriving 
eq. (|15|l from (JSJ. However, for a kernel given by a finite 
set of irreducible graphs, both BS QJ and LF (flT)l equa- 
tions are already approximate and it is not evident which 
of them is more " physical" . The physically transparent in- 
terpretation of the LF wave function makes it often more 
attractive. 



3 Ladder kernel 

We calculate here the kernel ^(7, z\ 7', z') of equation (|10|l 
for the ladder BS kernel, which reads: 



iK {L \k,k',p) 



i(-ig) 2 



(16) 



The singularities in the BS equation are removed by 
the analytical integration over (3. Equation (|10|l is valid We substitute it in eq. (|13|l . then substitute (|13|l in (|12|l 
for an arbitrary kernel iK, given by a Feynman graph. 



(k - fc') 2 - M 2 + «e 

I, then substitute 
and calculate the integrals. The details of these calcula- 
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tions are given in appendix |0 The result reads 

where W has the form: 

, , am 

W{j,z;j ,z ) 



< z' < z < 1 
-1 < z < z' < 1 
(17) 



(1-z) 2 



27r 7 + z 2 m 2 + (1 — z 2 )k 2 
(b+ - b-)(2b+b- - b+ - b-) 

(i-M(i-M 



+2b+b- log 



Ml -6+) 



(18) 



with a = 5 2 /(167rrn 2 ) and 



1 



6 = (1 - z)n 2 , 

bi = 7 + 7' — (1 — z)/i 2 — j'z — "fz 1 

+ (l-z') [z 2 m 2 + (1 - z V] , 
b 2 = -7(1 ~ z') 

- {z - z') [(1 - z)(l - z')k 2 + (z + z'- zz')m 2 ] . 

In the case /j, = (which constitutes the original Wick- 
Cutkosky model |3lTE2] ) we get, in particular, b = 6_ = 
and eq. H18|) obtains a more simple analytical expression 
that gives for the kernel: 



l/( L )( 7 ,z; 7 ',z') 



am 



1 



2tt 



7 + z 2 m 2 + (1 — z 2 )k 2 
1 



7' + z /2 m 2 + (1 - z /2 )k 2 
6{z-z') (l-z) 



7 + 7' hzr^h + -z 2 ^ 2 + (1 — z 2 )k 2 
6{z' ~ z) 



7 + 7' + z 2 to 2 + (1 - z 2 )k 2 




(1 + zO 



(19) 



We search for a solution of l|ll)|l in the form: 

g(l ) z) = 8( 7 )g(z). (20) 

The integration over 7' in both sides of equation (fTT)f> drops 
out. By setting 7' = everywhere, the kernel (|19fl turns 
into: 



L )(7,z; 7 ' = O,z') = 



1 



7 + z 2 m 2 + (1 — z 2 )k 2 
™2 



2tt 



z' 2 m 2 + (1 - z' 2 )n 2 



^4mz-z') + ^0(z'-z) 
(1 - z'j (1 + z') 



The prefactor 



7 + z 2 m 2 + (1 — z 2 )/c : 



is the same in both sides of equation IjlUfl and cancels. The 
7-dependence disappears thus from the equation which 
takes the simplified form: 



2tt 



dz' V(z,z')g{z') 



(21) 



with 



V(z,z') 



4 

(l-z') . 

(l±g) 

IT+FT' 



|(1 - z' 2 )M 2 

if -1 < z' < z < 1 
if -1< z < z' < 1 



(22) 



It exactly coincides with the Wick-Cutkosky equation 2 , 
3,12. 

Notice that in the [i — case the LF wave function 
(|14ll obtains the simple form |14U15j : 



x(l - x)g(l - 2.x) 



k\ + m 2 - x(l - x)M 2 



(23) 



with a known analytic dependence on k± variable. 



4 Numerical results 

Equation ijTUI) with the ladder kernel ((T7J) has been solved 
by using the same method as in |13j . i.e. by expanding the 
solution g on a spline basis |16j 



9(7, z) = ^2gijSt(l)Sj( z ) 



(24) 



over a compact integration domain fl — [0, "f m ax] x [— 1, +1] 
and validating the equation at some well chosen ensemble 
of collocation points {^i,Zj} C fl. The unknown coeffi- 
cient gij are determined by solving the resulting general- 
ized eigenvalue matrix equation 



A B{M)g = A(M)g 



(25) 



in which matrices B and A represent respectively the in- 
tegral operators of the left- and right-hand sides of I|1U|) . 
They both depend on the total mass of the system M 
through the parameter k defined in (|ll|l and the solution 
of the equation is provided by the values of M such that 
A(M) = 1. For the ladder kernel, the coupling constant a 
appears linearly in A and the problem can be formulated 
equivalently 

1 , . A(M) 
- B(M) g = g 
a a 
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e 

Fig. 1. Dependence of the coupling constant on the e param- 
eter of 1^ for At = 0.5 and B = 1.0. 

in which the inverse of the coupling constant appears as 
the eigenvalue of a linear system parametrised by M. 

It turns out that the discretized integral operator B 
has very small eigenvalues. They are unphysical but make 
unstable the solution of the system (|25() . To regularize B, 
we have added a small constant e to its diagonal part [T7| 
on the form: 

Bij fly + e Nij (26) 

where Nij is the equivalent of the Kronecker symbol 5ij 
in the bidimcnsional spline basis. This procedure allows 
us to obtain stable eigenvalues with an accuracy of the 
same order than e until values of e as small as 10~ 12 . 
We have plotted in fig. ^ the dependence of the coupling- 
constant on e for a system with binding energy B = 1.0 
and /i = 0.5. The convergence is very fast and a value 
e = 10~ 4 is enough to ensure a 4 digits stability on a. The 
real accuracy of a calculation is actually not determined by 
£ but rather by the grid parameters which were kept fixed 
in the results of fig-H These are essentially the value ^ ma x 
and the number of intervals N~ l and N z in each direction 
of Q. 



and /i = 0.5 the values displayed in tabled They corre- 
spond to "f m ax = 3, — 12, N z = 10. With all shown 
digits, they are in full agreement with the results we have 
obtained, similarly to jT^], by using the Wick rotation and 
the method of U|. Increasing e to 10~ 4 changes at most 
one unit in the last digit. This demonstrates the validity 
of our approach. 

We would like to remark the striking stability of the 
results with respect to iV 7 the number of grid points on 7. 
The value iV 7 = 12 used in our calculations was only for 
drawing purposes. In fact, the accuracy in calculating the 
eigenvalues in tabled is reached with N 7 = 1. This means 
that on the practical point of view our method leads to an 
equation whose solution is mostly one-dimensional and a 
number of grid points of the order of 10 on z-variable is 
enough to ensure an accuracy better than 10~ 4 . 

The weight function g for a system with /i = 0.5 and 
B = 1.0 is plotted in fig. [3 It has been obtained with 
£ = 10~ 4 and the same grid parameters than in table 
H Its 7-dependence is not monotonous and has a nodal 
structure; the ^-variation is also non trivial. We have re- 
marked a strong dependence of 5(7, z) relative to values 
of the £ parameter smaller than ~ 10 -4 , in contrast to 
high stability of corresponding eigenvalues. However the 
corresponding BS amplitude and LF wave function ip, 
obtained from 5(7,2) by the integrals © and (|14f> . show 
the same strong stability as the eigenvalues. 

The BS amplitude in Minkowski space in the rest frame 
p = is shown in fig. |3| The fc-dependence is rather 
smooth but the fco-dependence, due to poles of the prop- 
agators in ljT]l. exhibits a singular behaviour at 
fco = ± (fk ± if) i *- e - moving with fe and M. 

Note that our solution gives also the BS amplitude 
in Euclidean space, by substituting in @ &o = ik^. The 
Euclidean BS amplitude (E>E{ki,k) was found in this way 
in JH]. It is indistinguishable from the one obtained by a 
direct solution of the Wick-rotated BS equation. 

The corresponding LF wave function tp(k± , x) is shown 
in fig. 0] It is very similar to the LF wave functions dis- 
played in ref. |13j . though they obey different dynamical 
equations. It has a simpler structure than 3(7, z) in both 
arguments. 



Table 1. Coupling constant values as a function of the binding 
energy for ^ = 0.15 and fi = 0.5 obtained with y m ax = 3, 
iV 7 = 12, N z = 10 and e = 10" 6 . 



B 


a{fj, = 0.15) 


a(p = 0.50) 


0.01 


0.5716 


1.440 


0.10 


1.437 


2.498 


0.20 


2.100 


3.251 


0.50 


3.611 


4.901 


1.00 


5.315 


6.712 



By keeping e = 10 6 fixed and varying the grid param- 
eters to ensure four digits accuracy, we obtain for /i = 0.15 



5 Conclusion 

We have developed a method for solving the Bethe-Salpeter 
equation in the Minkowski space, i.e. without making use 
of the Wick rotation. 

The method is based on an integral transform of the 
original equation which removes its singularities. It is mo- 
tivated by the LF projection JJjJ of the BS amplitude. 

The transformed equation is formulated in terms of 
the weight function of the Nakanishi integral representa- 
tion 0, from which the original BS amplitude both in 
Minkowski and Euclidean spaces as well as corresponding 
LF wave function can be easily reconstructed. 

The equation has been obtained for scalar particles and 
applied to the ladder kernel. For zero-mass exchange the 



6 



V.A. Karmanov, J. Carbonell: Solving Bethe-Salpeter equation in Minkowski space 




y z 

Fig. 2. Function 17(7, 2) for fi — 0.5 and B — 1.0. On left - versus 7 for fixed values of z and on right - versus z for a few fixed 
values of 7. 




0.5 1 1.5 -1 -0.5 0.5 1 



k k Q 

Fig. 3. BS amplitude <&(ko, k) for /1 = 0.5 and B = 1.0. On left versus k for fixed values of ko and on right versus ko for a few 
fixed values of k. 

Wick-Cutkosky model is derived. For massive exchange, Acknowledgements 
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Fig. 4. Wave function ip(k±, x) for fi = 0.5 and B — 1.0. On left versus k± for fixed values of x and on right versus x for a few 
fixed values of k±. 



A Derivation of equation PU)! 

We substitute in equation the BS amplitude in the 
form © and apply to both sides the transformation ©. 
For the left-hand side this gives: 



dz' I di I dp (27) 
-i[l/A-{u-kf/{u-pf]g{i,z') 



7' + k 2 — k 2 — p ■ k z' — P' ( z' 



(ul-p) 



where (3' = (uj-p)(3. As seen from (|27|) . if) depends on three 
scalar pr 
relation 



scalar products k 2 , p ■ k and However, because of the 



p-k = 2 ( k 2 

ui ■ p 



1 



-M 



f2 



only two of them are independent. We use the LF variables 
J7|) and express through them the scalar products: 



k 2 ± 



4x{l - x) ' 



I, , f k 2 , +m 2 

p.k=-(l-2x)[^ r 

4 \ X \L — x ) 



M' 



oj ■ k 
uj ■ p 



By means of these relations, LF wave function (|27|l de- 
pends on k±, x. 

It is convenient to introduce other notations: 



7 



so that: 



1 - 2x, k 2 



(7 + z 2 m 2 ) 
1 -z 2 



-M 2 



p ■ k = 
uj ■ k 



z[y + z 2 m 2 + (1 — z 2 )k 2 



uj ■ p 



and 



4( 7 - 



2 ) k 2 +: 



1 



x(l — x) 



(28) 



(29) 



The integral (|27|) over j3' is simply calculated by means 
of the formula: 



dS 



-00 (J3x-y- ie) 3 



—S(x). 
y 



The result of transformation is given by eq. I|14|l . In terms 
of the variables 7, z it reads: 



ip(j,z) 



1 



(l-z 2 )g(y,zW 



7' + 7 + z 2 m 2 + k 2 (1 — z 2 ) 



( 3 °) 

Apart from the factor (1— z 2 )/(8y / 7r) (cancelled in the final 
equation) it is the left-hand side of eq. . Substitution 
of (J2J) and © in right-hand side of © results in the right- 
hand side of eq. (|10|l . 

The function ^(7,2) in eq. (|30|l is the usual two-body 
LF wave function. In terms of variables k±_ , x it obtains 
the form l|14|) . The normalization integral reads: 



No 



1 



(2tt) 3 
1 

(2^)3 



tp 2 (k±,x 
■ip 2 (j,z) 



d 2 k±dx 
' 2x{l - x) 
ird'ydz 
1-z 2 ' 



We would like to emphasize the mathematical nature 
of the above transformation. Namely, the integral trans- 
formation ©, applied to BS function J§J), obtains the form 
(|27|l . The function ip there still depends on the variables 
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k 2 ,p ■ k, like the initial BS function, but it is not singular. 
However, the variables 7', 2' in the integrand 5(7', z') and 
k 2 ,p ■ k in ip run over different domains: < 7' < 00, 
— 1 < z' < 1, whereas — 00 < k 2 < 0, —00 < p- k < 00. By 
eqs. (|28[) we replace the variables k 2 ,p- k by the new ones 
7, 2, which, by construction, vary in the same domain as 
7', z'. In these variables eq. I|27(l obtains the form (j20J), 
where the functions (7(7, z) and VKTj z ) before and after 
integration are now defined in the same domain, like it 
normally takes place in the integral equations. Therefore, 
in new variables 7,2 we will find equation for 5(7,2) in 
the domain of its definition. 

We separate the factor (1 — z 2 )/(8y / 7r), i.e. introduce 
ip related to tp as -0(7,2;) = (1 — z 2 )%j)(~f, z)/(8y / ir). That 
is: 

ni,Z)= 7- 72- (3i) 

Ja y + >f + z 2 m 2 + k 2 (1 - z 2 ) 

Eq. (|10|l can be rewritten for the function tp: 
1*00 1*1 

■tp(^f,z) = J dj' J dz iJ(7,z;7 / ,z')V'(7',^') I ( 32 ) 
where 

r-OO 

#(7, 2; 7', A = / <W( 7 , 2; 7", z')h{i\ 7', z') 



and hi^y" , 7', 2') is the kernel of the operator providing the 
relation inverse to eq. (|31|l . namely: 

/>oo 

Jo 

Numerical calculation of the kernel ^.(7, 7', z), then finding 
-ff(7, 2; 7', 2') and solving equation in the form (|32[l give 
the same results as for (|10() . 



B Derivation of the ladder LF equation (fT5)> 

We take, for a moment, uj = (1,0,0,-1), introduce the 
variables k'± = (k' x .k' y ), k'± = k' ± k' z and represent the 
integration volume d A k' in the right-hand side of eq. J5J) 
as 

d 4 k' = ^d 2 k' ± dk' + dk'_ = d 2 k' ± dx' {u-p)dj3', 

where x' is defined in J7J) and k'_ = oj^(3' (with w_ = 2). 
For arbitrary u> eq. (jHJ obtains the form: 



d 2 k' dx' 



(2tt)< 



J d(3${k + f3uj,p)= J d(3G 12 \k + f3<j,p) 
xiK(k + 0w, k' + (3'oj,p)<P{k' + I3'uj,p){oj ■ p)d0 (33) 

According to eq. BS amplitude ^(fc,p) contains as 
a factor the product of two free propagators J5J. We sep- 
arate the propagators, i.e., introduce the vertex function 

r(k, P y. 

<P(k,p) = G ( 12) (k 1 p)r(k,p) (34) 



and substitute 134JI in the right-hand side of Q33Jl ■ 

In order to derive LF equation (|15fl from (|33() . we 
should, integrating over /?', neglect the singularities of 
r(k' + (3'u>,p), i.e. deal with r vs. (3' as with a constant. 
That is, we will extract r from integral over /?' and in- 
troduce it back. This allows the following approximate 
transformation of the right-hand side of (|33|l (for short- 
ness we show only the /3, /^'-dependence and don't show 
integration over d 2 k' ± dx'): 

d(3G Q 12) (f3)K(0, p')G ( 12) ((3')r(f3')d(3' « 

r J dpG Q 12 \p)K{P>P')G 12 \(i'W = 

J G [p )dp J 



$G[ 12 \l3>W 

$dpG^\m{P,P)G < £ 2 \PW 



/G< 12J (W 



The integral J <P(l3')d(3' is understood as 
f^&ik' + P'uj,p)df3' and it is related by © to the LF 
wave function (as well as the left-hand side of 1)33(1 ). With 
explicit expression J2J for G Q 12 ^ we find: 

/>oo 

(oj-p) d(3G { 12) (k + (3Lu,p) 



x(l - x)(s — M 2 ) ' 



where s is defined in 129fl . With explicit expression (|16|l 
for K we obtain: 

poo 

(cu-p) 2 / d(3G ( Q 12) (k + (3u J ,p)K(k + f3Lu,k' + f3'uj,p) 



xG ( 12) (k' + p'uj,p)dp' 



(-4mVi£) 



x{l-x){s~M 2 ) V LF J x'(l- x')(s> - M 2 Y 

where V is the standard LF ladder kernel (see e.g. |10|1: 



V&\ k ' x ,x';k ± ,x,M 2 ) 
4:Tra8(x' - x) 



4:Tta9(x — x') 



and 



(x' ~ x){s a - M 2 ) (x — x')(sb — M 2 ) ' 



k\+ni 2 (k' j_ - k ± ) 2 + k'^+m 2 



(35) 



Sb 



x x' — x 1 — X 

k'l + m 2 (k' x - fc±) 2 + M 2 fei + m 2 



1-x 



In this way we derive the equation (|15fl with the kernel 
(H>J. 
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C Calculation of the kernel V^ L \j, z;j f , z'), 
eq. (HZD 

With the ladder kernel (|16fl the integral l|13|) obtains the 
form: 



I(k,p) 



-ilQirm a 

(27T) 4 



d 4 k' 



(k' - k) 2 - /.i 2 + ie 



1 



k' + p ■ k' z' — k 2 — 7' + ie 



where we put g 2 — l&irm 2 a. Using the formula: 



1 



3v 2 dv 



[a(l-v) + bvY 



and replacing k' by new integration variable q by the re- 
lation: 

k' = q + (1 — v)k — -vz'p, 



we get 
I(k,p) 



-ilQirm a 
(2^) 4 

2 /■! 



am 



3w 



v dv 



d 4 q 



[q 2 + A(p, k) + ie] 



with 



2n J [A(p,k) + te] 2 



A(p, k) = v(l — v)(k 2 + p ■ k z 1 ) — vm 2 
+ -v(l - vz' 2 )M 2 - (1 - u)/i 2 



(36) 



«7 



Now we make in (|36f) the replacement k — > k + /3a; 
and substitute it in l(T^|) . Scalar products k 2 ,p- k and also 
(lo ■ k)/(uj ■ p), which the kernel depends on, are expressed 
through the variables 7,2 by eqs. 1)28(1. Since 



A(p, k + {3uj) = A(p, k) + f3'v{\ - v)(z' - z) 
with (3 1 = (lo -p)f3, we obtain: 



at the value j3' ~ . . . + ie. We close the contour in the 
lower half-plane, i.e, take the residue at the pole of the 
first propagator: 







k 2 +p ■ k — k 2 



This gives: 



V {L \ 7 ,z; 1 ',z') = W(~ fl z; 1 ',z') 



with 

W{~(,z;i,z') 



am 2 (l — z) 2 



2tt 



7 + z 2 m 2 + (1 — z 2 )k 2 



1 „,2 



v dv 



and 



(38) 



D = v(l — v)(l - z')7 + v(l - z)i 
+ v(l - z)(l - z') 1 + z(l - v) + vz' 

(1-«)(1- z')z 2 +vz' 2 (l- z) 



+ (l-w)(l-z)/i 2 . 

In the case (it) the factor (z' — z) is positive and the 
pole in the variable j3' of the factor 1/ A + (3'v(l — v)(z' — 



z) +ie 



ie. We close the contour 



is at the value /3' ~ . 

in the upper half-plane, i.e, take the residue at the pole of 
the second propagator: 



k 2 — p ■ k — k 2 

P' = f- +«. 

l + z 



This gives: 



^( 7 ,z; T ',z') = %- Z ;7'-4 

with W defined in l|38l) . The integral for W is calculated 
analytically. In this way, we obtain eqs. (|1 7|) . I|18|l for the 
ladder kernel in the equation i|l(J|) . 



V^( 7 ,z- n ',z') = 



2tt 2 



< / d d» 
/o 



A(p, k) + (3'v(l - v)(z' - z) + ie 
1 



[k 2 + p ■ k - k 2 + (1 - z)(3' + ie] 
1 

[k 2 - p ■ k - k 2 - (1 + z)(3' + ie] ' 



(37) 



Both z and z' vary from —1 to 1. We consider two cases: 
(i) z' < z and (ii) z < z' . In the case (i) the factor [z' — z^j 
is negative and the (second order) pole in the variable [3' 



of the factor 1/ A + f3'v(l - v)(z' - z) + ie 



P7|) is 
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